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Abstract: Soil organic carbon (SOC) and its stable isotope composition reflect key information about the 
carbon cycle in ecosystems. Studies of carbon fractions in oasis continuous cotton-cropped fields can 
elucidate the SOC stability mechanism under the action of the human-land relationship during the 
oasification of arid land, which is critical for understanding the carbon dynamics of terrestrial ecosystems 
in arid lands under global climate change. In this study, we investigated the Alar Reclamation Area on the 
northern edge of the Tarim Basin, Xinjiang Uygur Autonomous Region of China, in 2020. In original 
desert and oasis farmlands with different reclamation years, including 6, 10, 18, and 30 a, and different soil 
depths (0—20, 20—40, 40—60 cm), we analyzed the variations in SOC, very liable carbon (Cyr), liable carbon 
(Cr), less liable carbon (Crx), and non-liable carbon (Cyr) using the method of spatial series. The 
differences in the stable carbon isotope ratio (6'C) and beta (f) values reflecting the organic carbon 
decomposition rate were also determined during oasification. Through redundancy analysis, we derived and 
discussed the relationships among SOC, carbon fractions, 8'C, and other soil physicochemical properties, 
such as the soil water content (SWC), bulk density (BD), pH, total salt (TS), total nitrogen (TN), available 
phosphorus (AP), and available potassium (AK). The results showed that there were significant differences 
in SOC and carbon fractions of oasis farmlands with different reclamation years, and the highest SOC was 
observed at the oasis farmland with 30-a reclamation year. Cvt, Cr, Cri, and Cyr showed significant 
changes among oasis farmlands with different reclamation years, and Cvr had the largest variation range 
(0.40-4.92 g/kg) and accounted for the largest proportion in the organic carbon pool. The proportion of 
Cyu in the organic carbon pool of the topsoil (0-20 cm) gradually increased. 8'3C varied from —25.61%0 to 
—22.58%0, with the topsoil showing the most positive value at the oasis farmland with 10-a reclamation year; 
while the f value was the lowest at the oasis farmland with 6-a reclamation year and then increased 
significantly. Based on the redundancy analysis results, the soil physicochemical properties, such as TN, AP, 
AK, and pH, were significantly correlated with Cr, and TN and AP were positively correlated with Cyr. 
However, 5'5C was not significantly influenced by soil physicochemical properties. Our analysis advances 
the understanding of SOC dynamics during oasification, revealing the risk of soil carbon loss and its 
contribution to terrestrial carbon accumulation in arid lands, which could be useful for the sustainable 
development of regional carbon resources and ecological protection in arid ecosystem. 
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1 Introduction 


The soil carbon pool is an important part of the terrestrial carbon pool, with approximately 1550 Gt 
C stored in the form of organic carbon (Lal, 2004). As a natural sink of atmospheric COs, soil 
organic carbon (SOC) pool plays an important role in coordinating soils, vegetation, and their 
relationships with the surrounding environment (Wiesmeier et al., 2019; Bossio et al., 2020). The 
accumulation and turnover rates of SOC significantly influence the global carbon budget and the 
stability of the climate system, thus impacting the global carbon cycle and climate change 
(Carvalhais et al., 2014). Land use is considered to be the most important factor affecting the 
accumulation and decomposition of SOC, and in arid and semiarid regions, the oasification process 
involving the gradual transformation of original desert to oasis farmlands has become a typical 
manifestation of land conversion efforts in these regions (Malik et al., 2018; Xue et al., 2019). In 
this process, agricultural activities, such as different crop cover types, tillage measures, and 
exogenous material inputs, can alter the plant rhizosphere processes and affect the stability of 
SOC, in turn changing the soil carbon storage (Gui et al., 2011; Finzi et al., 2015; Gong et al., 
2015). Therefore, obtaining a better understanding of SOC dynamics during oasification is 
essential for maintaining the stability of soil carbon pools in arid areas and the sustainable 
development of oasis farmlands. 

However, SOC, due to its high background value, shows a certain hysteresis in response to 
changes in land use methods and cannot optimally reflect changes in the SOC content or conversion 
rate over short periods of time (Post and Kwon, 2000). SOC has been found to be a heterogeneous 
mixture of compounds with different turnover cycles, and its chemical structure affects its stability 
(Kan et al., 2022). In general, SOC pool can be chemically divided into very liable carbon (Cyr), 
liable carbon (C1), less liable carbon (CL), and non-liable carbon (CNL), among which Cvr and Ci 
are more sensitive to land use changes than soil total organic carbon (Chan et al., 2001; Nath et al., 
2022). Several studies conducted in arid and semiarid regions have revealed that organic carbon 
fractions respond differently after land use changes, thus highlighting the unique sensitivity of these 
carbon fractions (Bremer et al., 1994; Yu et al., 2017; Zhong et al., 2019; He et al., 2021). For 
example, Cı and Cn in surface soils converted from original desert to croplands in the Yangi Basin 
of Xinjiang Uygur Autonomous Region, China, increased significantly (Zhang et al., 2014). 
However, little research has been conducted on oasification in arid and semiarid regions, and 
research on carbon fractions after the conversion of original desert into oasis farmlands with 
different reclamation years is especially scarce. The stability and dynamic changes that occur in soil 
carbon pool during the oasification process in arid regions are still unclear. 

Carbon stable isotope techniques can be used to describe the distribution of organic carbon, 
analyze the dynamics of SOC by studying the variation in stable carbon isotope ratio (6'3C), and 
interpret the turnover characteristics of organic carbon (Nissenbaum and Schallinger, 1974; 
Schaub and Alewell, 2009; Wang et al., 2018). In studies of 5!3C, the slope of the linear 
regression between the logarithm of SOC content and 8'°C is defined as beta (£) value; this term 
can be used as a representative of the organic carbon decomposition rate. The smaller the 6 value 
is, the faster the decomposition rate of SOC is (Zhao et al., 2019; Liu et al., 2021). During land 
use changes in arid or semiarid land areas, SOC content and distribution as well as 5'3C values 
vary with the degree of human disturbance (Liu et al., 2010; Sharafatmandrad, 2019; Zhu et al., 
2021). Therefore, it is important to investigate the trends and patterns of SOC and its fractions 
during the process of oasification to improve the soil carbon storage capacity, understand the 
stability of the soil carbon pool, and maintain the healthy development of oasis ecosystems. 

The oasis-desert ecotone at the northern edge of the Tarim Basin is the most representative oasis 
agricultural area and the main storage area of SOC in arid land. To support regional economic 
development, the exploitation of soil and water resources in this region has been expanding. The 
original desert has been gradually transformed into oasis farmlands through agricultural activities 
(e.g., tillage, irrigation, and fertilization), making the region an important area of continuous cotton 
cropping in China. Moreover, the stock of SOC has also changed in response to the conversion of 
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original desert to oasis farmlands (Zhang et al., 2003; Bai et al., 2021; Xue et al., 2022). Therefore, 
a deeper study of the stability of SOC pool during oasification is particularly important. 

To address this gap in knowledge, we selected original desert and oasis farmlands with 
different reclamation years to determine SOC and its fractions, 1C, and other soil 
physicochemical properties. The objectives of this study were to probe the changes in SOC and 
the distribution patterns of its fractions during oasification, investigate the changes in the °C 
and £ values that occur during oasification, and elucidate the effects of other physicochemical 
factors on SOC, carbon fractions, 5°C, and 2 values. Through this study, we aim to understand 
the ecological processes of soil carbon cycling during oasification in arid land to promote the 
sustainable development of regional carbon resources and provide a scientific basis for the study 
of the "carbon source and sink" effect in oases. 


2 Materials and methods 


2.1 Study area 


The Alar Reclamation Area is located at the northern edge of the Tarim Basin and the southern foot 
of the middle Tianshan Mountains, Xinjiang Uygur Autonomous Region of China. This region 
belongs to the alluvial plain area of the Tarim River, with slightly elevated terrain along the 
riverbanks and on both sides of the alluvial ditch; the terrain is also inclined from northwest to 
southeast. This area spans 281 km from east to west and 180 km from north to south, with an 
average elevation of 1012 m. Rainfall is scarce. Little snow falls in winter, and evaporation is 
strong on the land surface. The average annual precipitation ranges from 40.1 to 82.5 mm, and the 
average annual evaporation ranges from 1876.6 to 2558.9 mm. This area has good light and heat 
conditions and great temperature differences between day and night (Zhang, 2016). The local soil 
is alkaline and contains a high amount of salt (Peng et al, 2019). The main variety of planted cotton 
is "Xinluzhong 82", which is managed with the local integrated farming management methods, in 
which drip irrigation is applied 5—7 times during the reproductive period; the basal fertilizer 
applied comprises diammonium phosphate, urea, potassium sulfate, and mixed organic fertilizer. 


2.2 Sample collection 


Field soil samples were collected during June—July in 2020. In general, the closer the land area 
was to the inflow river, the earlier the oasis was formed. Therefore, old oases reclaimed as early 
as the 1950s and new oases reclaimed in recent years are distributed along the river to the inland 
region within the desert-oasis ecotone of the Alar Reclamation Area (Huang et al., 2021). We 
applied the spatial series method instead of the temporal series method to select four cotton fields 
with different reclamation years (i.e., 6, 10, 18, and 30 a) in a region with a similar geographical 
location and soil texture; irrigation was withheld from these fields for 15 d before sampling. The 
control group (CK) was an original desert with vegetation coverage consisting mainly of Alhagi 
sparsifolia, Tamarix spp., and Karelinia caspica (Ling et al., 2014). Each sample plot was 
approximately 0.12 hm?. After cotton harvest in October 2019, three typical sample squares with 
1 mx1 m areas were selected in each of the above sites with different reclamation years to collect 
soil samples by the quarter method. After removing obvious root residues, stones, and other 
debris, we crushed the soil samples and mixed them thoroughly before bringing them back to the 
laboratory for natural air-drying, grinding, and sieving to determine SOC content, carbon 
fractions, and other basic soil physicochemical properties. 


2.3 Laboratory analysis 


The content of SOC and its fractions was tested by the potassium dichromate method. According 
to the principle of the Walkley and Black method, the reaction heat produced by different 
concentrations of concentrated sulfuric acid resulted in different degrees of SOC oxidation 
(Walkley and Black, 1934; Chan et al., 2001). The carbon fraction oxidized by 6 mol/L sulfuric 
acid was regarded as Cyt. The difference in oxidizible organic carbon extracted between 9 and 6 
mol/L sulfuric acid were regarded as Cı. The difference in oxidizible organic carbon extracted 
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between 9 and 12 mol/L sulfuric acid were regarded as Cri. The carbon fraction oxidized by 12 
mol/L H2SO4 when compared with the total organic carbon was regarded as Cnt. 

Soil water content (SWC) was measured by the drying method. The soil samples were dried to 
a constant weight at 105°C, and the wet and dry weights of the soil samples were recorded. We 
used the cutting ring method to determine bulk density (BD). The original-state soil samples taken 
by the cutting ring were brought back to the laboratory and dried to a constant weight at 105°C. 
After cooling, the dry weights were measured, and BD were calculated (Huntington et al., 1989). 
A total of 20-g soil was dried and passed through a 1-mm sieve, and 100-mL distilled water was 
added for shock filtration. Soil pH was determined by using a pH meter (PHS-3C, Shanghai 
Yidian Scientific Instrument Co. Ltd., Shanghai, China). Total salt (TS) content was the sum of 
the eight ions' concentrations (including C037, HCO3,, Cl, S047, Ca?*, Mg?*, Nat, and K*). 
CO3* and HCO; were determined by the double indicator neutral method; CI- was determined 
by the AgNO; (Mohrs's titration) method; SO4% was determined by indirect Ethylene Diamine 
Tetraacetic Acid (EDTA) titration; Ca?* and Mg** were determined by complexometric titration 
with EDTA; and Na* and K* were determined by flame photometry. Total nitrogen (TN) was 
determined by the Kjeldahl method (K370, BUCHI Labortechnik AG, Flawil, Switzerland). 
Available phosphorus (AP) was extracted with NaHCO; and determined by colorimetry. Available 
potassium (AK) was determined by the flame photometry method. 

The expression of ô!?C is as follows (Breecker et al., 2015): 


13 12 
soa ECC) same 1), 1000% (1) 
rer 7 0) eee 


where 8!3C is stable carbon isotope ratio (%o), which is used to describe the variation in the °C 
when a sample is compared to a standard sample (the error of 5!3C is +0.2%o); (3C/!*C) sample is the 
ratio of '3C to °C in the sample; and (13C/!?C)standara is the ratio of '3C to !?C in standard material 
of Pee Dee Belemnite (PDB) in South Carolina, USA, which is defined as '3C/!*C=0.01124 
(O'Leary, 1981). 


2.4 Statistical analysis 


The data obtained in this study were first organized and analyzed in Excel 2016 for the basic 
statistical analysis; then, analysis of variance and least significant difference were used to test 
SOC, carbon fractions, 5'°C, and other physicochemical factors in different reclamation years and 
soil depths and to analyze the interaction between the soil depth and reclamation years using 
SPSS 22.0 software. The Kruskal-Wallis test was used to compare the variability in J values 
among different reclamation years. Canoco 5.0 software was selected for the redundancy analysis 
(RDA) to investigate the correlations among carbon fractions, 5'°C, and other physicochemical 
factors of the soil. Pearson correlation analysis was used to analyze the correlations among SOC, 
513C, and soil environmental factors. Origin 2020 and R Studio were used for plotting. 


3 Results 


3.1 Variations of soil physicochemical properties during oasification 


The physicochemical properties of the soil samples are presented in Table 1. SWC measured at 
different soil depths remained the lowest in original desert (CK), with BD ranging from 1.25 to 
1.43 g/cm?. The value of pH was the highest in original desert (8.66). The soil was alkaline, with 
a decreasing trend observed alongside continuous tillage. Variations in TS was large, ranging 
from 1.59 to 18.31 g/kg. TN, AP, and AK contents tended to be highest at the soil depth of 0-20 
cm after 6-a reclamation. The effects of reclamation year, soil depth, and their interaction on soil 
physicochemical properties are shown in Table 2. Soil physicochemical properties except AK 
were significantly affected by reclamation year. Soil depth had significant effects on SWC, BD, 
TS, TN, and AP. Among them, BD, TN, and AP were significantly affected by the interaction 
between reclamation year and soil depth (P<0.05). 
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Table 1 Variations of physicochemical properties in different reclamation years and soil depths 


Reclama- Soil depth 


; SWC (%) BD (g/cm*) pH TS (g/kg) TN (g/kg) AP (mg/kg) AK (mg/kg) 
tion year (cm) 

0-20 10.62+1.50 1.37+0.04 8.6640.14  18.31414.26 0.25+0.03 2.9940.35 164.53419.42 

CK 20-40 12.18+1.60 1.36+0.05 8.39+0.14 9.9841.65  0.2040.04  3.7340.95 170.22+18.46 


40-60 13.6941.73 1.4340.05 8.31+40.04 6.28+0.73  0.24+0.03 3.59+40.50 162.49422.87 
0-20 18.68+1.32 1.29+0.04 8.0740.05 6.5343.54 0.37+0.03  18.3345.89 185.81+31.15 
6-a 20-40 24.96+1.79 1.27+40.03 8.04+0.03 2.2341.26 0.32+0.03 19.8544.94 177.88+16.29 
40-60 26.25+1.41 1.40+40.04 8.1540.11 1.59+0.50 0.26+0.03 11.15+2.09 171.90+421.23 
0-20 21.2941.57 1.25+40.02 7.96+0.12 7.41+4.68 0.56+0.03 28.2743.99 229.31+49.75 
10-a 20—40 26.09+1.49 1.2740.03 7.9840.03 2.5140.71  0.4640.07 21.2044.45 201.76+96.14 
40-60 28.69+2.61 1.28+40.04 8.01+40.07 1.9840.80  0.37+0.09 17.84+3.95 211.91+65.23 
0-20 20.8840.99 1.28+0.02 7.92+0.21 14.0145.37  0.55+0.06 30.5345.09 213.26+23.62 
18-a 20-40 25.4941.53 1.27+40.05 7.82+0.23 6.76+2.77 0.5440.04 22.80+5.41 227.52+27.83 
40—60 28.59+2.35 1.26+0.02 7.8140.19 2.14+0.22 0.4140.05  14.8542.11 211.34+17.06 
0-20 21.8041.56 1.32+0.04 7.80+0.13 16.29+5.33  0.58+0.05 31.30+1.66 211.08+42.75 
30-a 20-40 24.62+1.15 1.3340.02 7.81+0.04 8.3543.85 0.4740.09 21.54+4.17 208.82+31.43 
40-60 24.72+1.69 1.36+0.02 7.82+0.17 4.24+0.53  0.33+0.04 15.51+2.80 205.05+20.43 


Note: CK, original desert; SWC, soil water content; BD, bulk density; TS, total salt; TN, total nitrogen; AP, available 
phosphorus; AK, available potassium. Mean+SD. 


Table 2 Effects of reclamation year, soil depth, and their interaction on soil physicochemical properties 


Source Controlled variable Sum of square df Mean square F value P value 
Reclamation year 1093.22 4 273.31 98.32 <0.01 
SWC Soil depth 259.77 2 129.88 46.72 <0.01 
Reclamation yearxSoil depth 44.19 8 5.52 1.99 0.08 
Reclamation year 0.07 4 0.02 14.35 <0.01 
BD Soil depth 0.01 2 0.01 4.68 0.02 
Reclamation yearxSoil depth 0.02 8 0.01 2.56 0.03 
Reclamation year 2.40 4 0.60 35.73 <0.01 
pH Soil depth 0.05 2 0.02 1.37 0.27 
Reclamation yearx Soil depth 0.20 8 0.03 1.51 0.19 
Reclamation year 443.60 4 110.90 5.13 0.03 
TS Soil depth 679.98 2 339.99 15.74 <0.01 
Reclamation yearx Soil depth 84.58 8 10.57 0.49 0.85 
Reclamation year 0.49 4 0.12 49.93 <0.01 
TN Soil depth 0.14 2 0.07 29.75 <0.01 
Reclamation yearx Soil depth 0.06 8 0.01 3.18 0.01 
Reclamation year 2506.30 4 626.58 46.14 <0.01 
AP Soil depth 705.92 2 352.96 25.99 <0.01 
Reclamation yearx Soil depth 343.52 8 42.94 3.16 0.01 
Reclamation year 19450.03 4 4862.51 3.07 0.31 
AK Soil depth 514.99 2 257.50 0.16 0.85 
Reclamation yearx Soil depth 1563.37 8 195.42 0.12 0.99 


Note: df is degree of freedom. 
3.2 Variations of soil organic carbon (SOC) and carbon fractions during oasification 


SOC varied significantly during the process of oasification (Table 3). At different soil depths, 
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SOC was the highest at the oasis farmland with 10-a reclamation year. Cv, had a wide range of 
variation (0.40-4.92 g/kg) and showed the highest accumulation at the soil depth of 20-40 cm at 
the oasis farmland with 30-a reclamation year. The highest Cnr (3.62 g/kg) was also observed at 
the oasis farmland with 30-a reclamation year, and the lowest was observed in original desert 
(CK). Cv: changed significantly among different reclamation years and soil depths, while Cu, 
Cit, and Cyne exhibited significant changes only among different reclamation years (Table 4). 
SOC was significantly correlated with Cv (R7=0.94; P<0.01), CL (R?=0.42; P<0.01), Ci 
(R?=0.38; P<0.01), and Cn (R?=0.88; P<0.01) in all reclamation years and soil depths (Fig. 1). At 
the soil depths of 0-20 and 20—40 cm, SOC pool were dominated by liable fractions (i.e., Cvt, CL, 
and Cı) (Fig. 2). At the soil depth of 0-20 cm, Cni (18.63%-—34.52%) showed a gradual increase 
with increasing reclamation years. At the soil depth of 40-60 cm, Cnr accounted for the main 
proportion (56.30%) of each carbon fraction in original desert. 


Table 3 Variations of soil organic carbon (SOC) and carbon fractions in different reclamation years and soil depths 


Reclamation year Soil depth (cm) Cvi (g/kg) C (g/kg) Cit (g/kg) Cnt (g/kg) SOC (g/kg) 
0-20 0.66+0.12 0.21+0.05 0.40+0.08 0.29+0.07 1.44+0.21 

CK 20-40 0.53+0.17 0.19+0.01 0.24+0.14 0.36+0.18 0.93+40.24 
40-60 0.40+0.08 0.08+0.01 0.12+0.03 0.77+0.33 1.20+0.32 

0-20 1.56+0.06 0.72+0.30 0.08+0.02 0.60+0.28 2.99+0.20 

6-a 20—40 1.04+0.14 0.28+0.08 0.04+0.01 0.64+0.01 1.88+0.38 
40-60 0.61+0.05 0.16+0.08 0.05+0.02 0.48+0.07 1.3040.18 

0-20 2.13+40.39 0.56+0.36 0.45+0.20 0.72+0.28 3.86+0.83 

10-a 20-40 1.46+0.60 0.74+0.44 0.19+40.09 0.5140.18 2.90£1.12 
40-60 0.82+0.05 0.64+0.24 0.32+0.18 0.29+0.07 2.07+0.08 

0-20 1.80+£0.28 0.84+0.28 0.13+0.05 0.96+0.14 3.67+0.21 

18-a 20-40 1.24+1.07 1.44+0.90 0.29+0.18 0.66+0.28 3.24+1.02 
40-60 0.74+0.70 0.59+£0.33 0.3740.14 0.37+0.12 2.074£1.75 

0-20 4.34+0.47 1.97+0.41 0.56+0.08 3.62+0.59 10.48+0.59 

30-a 20-40 4.92+0.64 0.9840.51 0.64+0.14 3.08+0.53 9.63+0.53 


40-60 0.60+1.12 1.49+0.48 0.64+0.14 2.15+1.19 7.7442.15 


Note: Cyr, very liable carbon; Cy, liable carbon; Cy, less liable carbon; Cwi, non-liable carbon. Mean+SD. 


Table 4 Effects of reclamation year, soil depth, and their interaction on carbon fractions 


Source Controlled variable Sum of square df Mean square F value P value 
Reclamation year 76.03 4 19.01 79.73 <0.01 
Cyt Soil depth 6.27 2 3.13 13.14 <0.01 
Reclamation yearxSoil depth 2.64 8 0.33 1.38 0.24 
Reclamation year 9.57 4 2.39 21.07 <0.01 
CL Soil depth 0.54 2 0.27 2.38 0.11 
Reclamation yearx Soil depth 2.66 8 0.33 2.93 0.02 
Reclamation year 1.44 4 0.36 8.50 <0.01 
Cit Soil depth 0.02 2 0.01 0.18 0.84 
Reclamation yearxSoil depth 0.31 8 0.04 0.93 0.51 
Reclamation year 41.56 4 10.39 47.35 <0.01 
Cnr Soil depth 1.35 2 0.68 3.08 0.06 
Reclamation yearxSoil depth 3.18 8 0.40 1.81 0.12 
Reclamation year 370.84 4 92.71 118.04 <0.01 
SOC Soil depth 19.44 2 9.72 12.38 <0.01 
Reclamation yearxSoil depth 6.03 8 0.75 0.96 0.48 
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Fig. 1 Relationships of soil organic carbon (SOC) with different carbon fractions. Cv, very liable carbon; CL, 
liable carbon; CLL, less liable carbon; Cwi, non-liable carbon. 
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Fig. 2 Proportions of Cv, CL, Cit, and Cnr in total SOC at the soil depth of 0-20 cm (a), 20—40 cm (b), and 
40—60 cm (c). CK, original desert. 
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3.3 Variations of stable carbon isotope ratio (6'°C) and beta (£) value during oasification 


As shown in Figure 3, 6!3C varied significantly (from —25.61%o to —22.58%o) with reclamation 
years and soil depths. 5'°C value at the soil depths of 0-20 and 40-60 cm changed significantly 
(P<0.05) among different reclamation years. Among them, at the soil depth of 0-20 cm, &!3C 
showed a significant increase at the oasis farmland with 10-a reclamation year. At the soil depth 
of 20-40 cm, the variation in 6!3C in each reclamation period fluctuated but was not significant. 
At the soil depth of 40—60 cm, 5'°C showed significant decreases at the oasis farmlands with 6-a 
and 10-a reclamation year and significant increases at the oasis farmlands with 18-a and 30-a 
reclamation year (P<0.05). 5°C changed significantly among different soil depths at the oasis 
farmland with 10-a reclamation year (P<0.05). In addition, £ values ranged throughout the 
process of oasification in oasis farmlands with the following rank: 18-a>10-a>30-a>CK>6-a. In 
addition, there was a significant increase from the oasis farmlands with 6-a reclamation year to 
the oasis farmlands with 10-a reclamation year (Fig. 4). 
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Fig. 3 Changes of stable carbon isotope ratio (8!°C) in different reclamation years and at different soil depths. 

Different lowercase letters represent the significant differences among the three soil depth in the same 

reclamation year (P<0.05), and different uppercase letters represent the significant difference among different 
reclamation years at the same soil depth (P<0.05). Bars mean standard errors. 
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Fig. 4 Changes of beta (£) value in different reclamation years. * denotes significant at P<0.05 level based on 
the Kruskal-Wallis (KW) test. The top, middle, and bottom lines of the box represent the upper quartile, median, 
and lower quartile, respectively; the black dot in the box represents the mean; and the top bar and bottom bar 
represent the maximum and minimum, respectively. 
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3.4 Effect of soil physicochemical properties during ossification 


The ordination relationships among SOC, carbon fractions, 5!°C, and other soil physicochemical 
properties are shown in Figure 5. Among them, TN, pH, AP, and AK exhibited the longest 
connecting lines, and the cosine value of the angle between Cy, and Cr was relatively large. These 
results suggested that TN, pH, AP, and AK have large correlation coefficients with carbon fractions 
and are positively correlated, except for pH. TN, AP, AK, and pH greatly influenced on the 
transformation of SOC. The connecting lines of SWC and TS were long, but the cosine value of 
this angle was relatively low. SWC and TS were positively correlated with carbon fractions, but 
these correlations were not obvious. Low cosine values were observed between 8'°C and soil 
physicochemical properties, though these correlations were not significant. The effects of seven 
soil physicochemical properties on carbon fractions differed. The results showed that the order of 
importance was TN>pH>AP>AK>TS>SWC>BD (Table 5). Among them, the effects of TN, AP, 
pH, and AK on carbon fractions reached significant levels (P<0.05), while the effects of TS, SWC, 
and BD on SOC, carbon fractions, and 8'C did not reach the significant level (P>0.05; Table 6). 
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Fig. 5 Redundancy analysis (RDA) of SOC, carbon fractions, and 8'°C with soil physicochemical properties. 
SWC, soil water content; BD, bulk density; TS, total salt; TN, total nitrogen; AP, available phosphorus; AK, 
available potassium. 


Table 5 Importance and significance level of soil environmental factors to SOC, carbon fractions, and 6!3C 


Soil environmental factor Rank Proportion (%) F value P value 
TN 1 26.1 4.6 0.016 

pH 2 25.5 4.5 0.024 

AP 3 24.3 4.2 0.020 

AK 4 20.4 3.3 0.036 

TS 5 12.1 1.8 0.162 

SWC 6 6.4 0.9 0.432 

BD 7 3:1: 0.5 0.666 


4 Discussion 


4.1 Effect of oasification on SOC 
In this study, different reclamation years were found to significantly affect SOC, and a significant 
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Table 6 Correlation analysis of SOC, carbon fractions, and 5!3C with soil physicochemical properties 


SWC BD pH TS TN AP AK 

SOC 0.86 —0.51 —0.76 0.29 0.96 0.93 0.45 

Cv 0.16 —0.19 —0.49 0.38 0.66” 0.63" 0.44 
CL 0.37 —0.28 -0.72™ 0.23 0.70 0.66" 0.68" 

Cit 0.11 0.00 —0.36 0.30 0.35 0.22 0.48 

Cnt 0.13 0.14 -0.51 0.35 0.44 0.45 0.31 

8BC 0.38 —0.29 —0.09 —0.47 0.09 0.25 0.12 


Note: * and “ denote significant at P<0.05 and P<0.01 levels, respectively. 


increase in SOC was observed in the whole soil profile after reclamation from the original desert 
to continuously cropped farmlands (Table 3). This finding is consistent with the conclusions of 
previous research that the carbon sequestration process that occurs in arable lands under 
agricultural tillage is ofa greater magnitude than that which occurs in uncultivated soils (DeGryze 
et al., 2004). The reason for this difference may be that crop residues are effectively accumulated 
in soils during the reclamation and cultivation of oasis farmlands (Rakkar and Blanco-Canqui, 
2018). At the same time, the large amount of organic fertilizer applied impacts the physical and 
chemical processes of the soil and accelerates biological processes within the soil, thereby 
accelerating the biodegradation of organic carbon (Vázquez et al., 2016; Zhou et al., 2022). Crops 
absorb water and nutrients from the soil through their root systems while growing and 
developing, so the humus layer of soil in oasis farmlands is thicker than that in the original desert 
(Jobbágy and Jackson, 2000; Nakamoto et al., 2012; Zhang et al., 2014). The increased root 
biomass promotes the formation and stability of soil aggregates, and the increased stability of the 
soil structure contributes to the accumulation of SOC (Li et al., 2020; Okolo et al., 2020). 
In this study, SOC and Cyr were found to be highly and significantly influenced by soil depth. 
With the increase of soil depth, SOC and Cr decreased (Table 3). When a large amount of litter 
accumulates and is converted into organic matter on the soil surface, the carbon source 
supplementation obtained in deeper soil layers is reduced, and the microbial activity is also 
reduced, thus affecting the contents of more active carbon fractions, such as Cyr, Ci, and Cız; 
this process is consistent with the findings of many other studies performed in arid areas (Cotrufo 
et al., 2015; Joly et al., 2020). The proportions of these more active carbon fractions in the total 
organic carbon reflect the quality and stability of SOC to a certain extent. In the results, the 
proportions of carbon fractions, except Cn, in the topsoil decreased during the process of 
oasification (Fig. 2), indicating that the deepening of reclamation could promote the conversion of 
active organic carbon to passive organic carbon and thereby improve the stability of carbon pool. 
In the SOC pool of oasis farmlands, SOC was significantly correlated with Cvi, CL, and Ci 
(Fig. 1). The proportion of Cnr at the soil depth of 0-40 cm was lower than the proportion of the 
sum of the active carbon fractions (Cv, Ci, and Cir) (Fig. 2). This finding was not consistent 
with the results of previous studies (Rovira et al., 2012; Zhang et al., 2014; Ahmed et al., 2016). 
Original desert in arid lands has low initial SOC. When original desert is reclaimed into oasis 
farmlands, the input of organic matter provides a sufficient carbon source for microbial growth 
and reproduction and accelerates the decomposition of organic matter in soil surface (Prommer et 
al., 2019). Because the sources of soil active organic carbon are microorganisms, which are active 
in crop residues, the increased microbial metabolic functions and secretions are important reasons 
for the increased carbon fractions being more easily decomposed during the conversion of 
original desert into oasis farmlands (Rovira and Vallejo, 2002). Compared to carbon, which is 
difficult to decompose, active organic carbon has a higher turnover rate and is more sensitive to 
soil physicochemical properties (e.g., pH and SWC), suggesting that the oasification process by 
which original desert is conversion to oasis farmlands increases SOC but also risks carbon loss 
(Benbi et al., 2015; Guan et al., 2022). It is noteworthy that in surface soils, although Cnr 
accounts for a relatively low proportion of the total SOC (Fig. 2), this fraction has a gradually 
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increasing trend, revealing the promotional effect of the oasification process on the terrestrial 
carbon storage capacity in arid lands. 


4.2 Effect of oasification on soil 6°C and £ value 


The organic carbon isotope composition is extremely weakly fractionated during the 
decomposition of plant residues and the layer-by-layer accumulation process (Wang et al., 2017; 
Li et al., 2022). The isotope values obtained in this study ranged from —25.61%o to —22.58%o (Fig. 
3), slightly higher than the average isotope value of C3 plants; our result was similar to the 
findings of other studies (de Rouw et al., 2015; Chen et al., 2018), indicating that 5'3C becomes 
enriched in the studied soils. The reason for this may be the fractionation of organic carbon 
isotopes. Compared to C, '7C more easily enters the CO2 produced by the decomposition of soil 
organic matter, and '°C increases 8'°C value in the soil by forming SOC (Li and Schaeffer, 2020). 
At different soil depths, &!3C values at the oasis farmlands with 6-a and 10-a reclamation year 
were significantly positive compared to other oasis farmlands (Fig. 3); this finding may have been 
related to the improvements in microbial activity and fecundity caused by soil environmental 
changes that occur during the tillage process, as these improvements are conducive to the 
degradation of refractory substances in crops (Wynn et al., 2006; Brunn et al., 2014; Xiao et al., 
2019; Sokol et al., 2022). However, 85°C showed a significant negative tendency at the oasis 
farmlands with 18-a and 30-a reclamation year (Fig. 3), indicating that the degradation of 
microorganisms has a certain limit. The reason for this may be that agricultural activities such as 
tillage and irrigation can accelerate infiltration and precipitation (Blaser and Conrad, 2016). With 
the input of crop litter to soil surface during the tillage process, lignin, cellulose, and other 
refractory substances accumulate in the soil, resulting in a decrease in 8'C (Gao et al., 2016). 
Meanwhile, £ value can indicate the decomposition rate of SOC (Zhang et al., 2021). In the 
results, 6 value at the oasis farmland with 6-a reclamation year was the lowest and changed 
significantly compared with that at the oasis farmland with 10-a reclamation year (Fig. 4), 
indicating that the main reason for the positive 5'°C values measured at the oasis farmland with 
6-a reclamation year was the increase in the decomposition rate. However, the decomposition rate 
of the samples at the oasis farmland with 10-a reclamation year was low but was still positive, 
which may have been related to the high SWC at the oasis farmland with 10-a reclamation year 
accelerating the migration and accumulation of organic carbon in the soil. 


4.3 Interaction between soil carbon and soil physicochemical properties 


Microorganisms secrete extracellular enzymes that participate in the metabolic processes of soil 
organic matter, thus procuring energy for themselves. When the demand for carbon sources 
increases, the related enzyme activities increase (Trivedi et al., 2016; Calabrese et al., 2022). In 
this study, Cv and Cy had significant correlations with TN (Fig. 5); this may have been related to 
the application of fertilizers during tillage. On the one hand, nitrogen in soil contributes to 
condensation with lignin and other substances in crop residues, inhibits the activity of enzymes in 
decomposing passive organic carbon, and forms heterocyclic substances (such as indole) and 
phenols, which are difficult to decompose, thus promoting the accumulation of soil passive 
organic carbon (Ma et al., 2022). On the other hand, nitrogen can increase the production of 
organic matter such as lignin and other materials that are difficult to degrade, which may cause 
elevated soil passive organic carbon contents and distribution ratios, thus improving the stability 
of the soil carbon pool (Chen et al., 2018). AP and AK in soil samples showed significant 
increases throughout the process of oasification and were highly significantly correlated with Cvt 
or Ci (Fig. 5; Table 5), indicating that an increase in phosphorus in soil may have a stimulating 
effect leading to the mineralization of organic carbon. Therefore, the nutritional status of oasis 
farmlands can be improved, and this improvement is conducive to the absorption and utilization 
of soil nutrients by crop roots (Nottingham et al., 2015; Tan et al., 2017; Woittiez et al., 2019; Lu 
et al., 2020). 

In arid areas, the salt in groundwater tends to move upward to the crop rhizosphere with the 
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evaporation of soil water and easily accumulates in the soil surface, resulting in alkaline soils 
(Zhang et al., 2014). The soil pH affects the mineralization and decomposition of SOC mainly by 
influencing microbial and enzyme activities and therefore has an inhibitory effect at high pH 
values (pH>8.5), decelerating the turnover of SOC (Malik et al., 2018). With the reclamation of 
oasis farmlands, a significant decrease in the soil pH (from 8.66 to 7.80) and a significant 
negative correlation between the soil pH and Cı can be observed (Tables 1, 2, and 6; Fig. 5), 
suggesting that a decline in soil alkalinity contributes to increases in microbial activities and the 
reproduction rate, which in turn benefit mineral adsorption by microbial residues and improve the 
decomposition of organic matter, resulting in the increased proportion of the active carbon pool 
throughout the oasification process. 


5 Conclusions 


SOC and its stable isotope composition were studied by the method of spatial series instead of 
temporal series in soil samples with different reclamation years in Alar Reclamation Area, 
Xinjiang Uygur Autonomous Region of China. The effects of soil physicochemical properties 
were analyzed to better understand the soil carbon dynamics during oasification process. The 
results showed that SOC and carbon fractions of oasis farmlands varied significantly among the 
regions with different reclamation years, and the highest SOC was observed at the oasis farmland 
with 30-a reclamation year. With the increase of soil depth, SOC reduced. Cvr had the largest 
variation range (0.40—4.92 g/kg) and accounted for the largest proportion in the SOC pool. The 
proportion of carbon fractions except Cnr decreased during oasification, but still occupied the 
main component. However, the proportion of Cnt in the carbon pool of topsoil (0-20 cm) 
gradually increased. 5'°C varied from —25.61%o to —22.58%o, with the topsoil showing the most 
positive value at the oasis farmland with 10-a reclamation year; the 2 value was the lowest at the 
oasis farmland with 6-a reclamation year and then increased significantly. The redundancy 
analysis results showed that TN and AP were positively correlated with Cyt. The main soil 
physicochemical properties influencing C were TN, AP, AK, and pH. Oasification process, the 
conversion from original desert to oasis farmlands, altered the physicochemical properties of 
natural arid land soil, contributing to the storage of organic carbon in the arid ecosystem in 
northwestern China. These results provide a scientific basis for the study of the "carbon source 
and sink" effect in oases. 
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